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Abstract 

We consider rate swaps which pay a fixed rate against a floating rate in 
presence of bid-ask spread costs. Even for simple models of bid-ask spread 
costs, there is no explicit strategy optimizing an expected function of the 
hedging error. We here propose an efficient algorithm based on the stochas¬ 
tic gradient method to compute an approximate optimal strategy without 
solving a stochastic control problem. We validate our algorithm by numerical 
experiments. We also develop several variants of the algorithm and discuss 
their performances in terms of the numerical parameters and the liquidity 
cost. 


1 Introduction 

Classical models in financial mathematics usually assume that markets are per¬ 
fectly liquid. In particular, each trader can buy or sell the amount of assets he/she 
needs at the same price (the “market price”), and the trader’s decisions do not 
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affect the price of the asset. In practice, the assumption of perfect liquidity is 
never satished but the error due to illiquidity is generally negligible with respect 
to other sources of error such as model error or calibration error, etc. 

However, the perfect liquidity assumption does not hold true in practice for in¬ 
terest rate derivatives market; the liquidity costs to hedge interest rate derivatives 
are highly time varying. Even though there exist maturities for which zero-coupon 
bonds are liquid, bonds at intermediate maturities may be extremely illiquid (no¬ 
tice that the underlying interest rate is not directly exchangeable). Therefore, 
hedging such derivatives absolutely needs to take liquidity risk into account. In this 
context, dehning and computing efficient approximate perfect hedging strategies 
is a complex problem. The main purpose of this paper is to show that stochastic 
optimization methods are powerful tools to treat it without solving a necessarily 
high dimensional stochastic control problem, under the constraints that practi¬ 
tioners need to trade at prescribed dates and that relevant strategies depend on a 
hnite number of parameters. More precisely, we construct and analyze an efficient 
original numerical method which provides practical strategies facing liquidity costs 
and minimizing hedging errors. 

The outline of the paper is as follows. Section 2 introduces the model. In 
Section 3, we present our numerical method and analyze it from a theoretical 
point of view within the framework of a Gaussian yield curve model. Section 4 
is devoted to a numerical validation in the idealistic perfect liquidity context. In 
Section 5, we develop an empirical study of the efficiency of our algorithm in the 
presence of liquidity costs. 


2 Our settings: swaps with liquidity cost 

2.1 A short reminder on swaps and swaptions hedging with¬ 
out liquidity cost 

One of the most common swaps on the interests rate market is as follows. The 
counterparts exchange two coupons; the hrst one is generated by a bond (with 
a constant hxed interest rate) and the second one is generated by a floating rate 
(e.g. a LIBOR). 

Definition 1. In a perfectly liquid market, the price at time t of a zero-coupon 
bond paying 1 at time T is denoted hy B{t, T). The linear forward rate L{Tf, Tb, Te) 
is the fair rate decided at time Tp determining the amount at time Te obtained 


2 



by investing 1 at time Tb- The following relation is satisfied: 


L{Tb,Tb,Te) 


1 ( B{Tb,Tb) \ 
Te-Tb\B{Tb,Te) )' 


( 1 ) 


A swap contract specifies: 


• an agreement date t 

• a time line (t <) Tq < ■ ■ ■ < T^ 

• a fixed interest rate r 

• a floating interest rate 

• the payoff at each time Tj (1 < f < N), that is, 

P{t) := (T, - T,_i) {r-L (T,_i, T,_i, T,)). 
From (1), we deduce the equivalent expression 


P{z) =r{T,- T,_i) 


1 


+ 1 . 


( 2 ) 

(3) 


In the sequel, we consider that the hxed rate r is chosen at the money (thus the 
swap at time t has zero value), and that the swap hxed coupons are received by 
the trader. 

In the idealistic framework of a market without liquidity cost, the trader buys 
or sells quantities of zero-coupon bonds at the same price (i.e. the market price), 
and there exists a discrete time perfect hedging strategy which is independent of 
any model of interest rates. In view of (3), the replication of the payoff P{i) at 
time Ti can be split into three parts: 


• the hxed part r{Ti — Tj_i) is replicated statically at time t by selling r{Ti — 
Tj_i) zero-coupon bonds with maturity Tj. 

• the boating part l/i?(Tj_i, T*) is replicated dynamically at time Ti_i by 
buying l/B{Ti_i,Ti) zero-coupon bonds with maturity Tj. The price of this 
transaction is equal to 1. 

• the last (hxed) part 1 is used at time T* to buy l/B(Ti,Ti+i) zero-coupon 
bonds with maturity Tj+i. 
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It is easy to see that this strategy is self-hnancing at times Ti, • ■ ■ , Tjv_i. To make 
it self-hnancing at time t also, at this date one buys 1 zero-coupon bond with 
maturity Tq and sells 1 zero-coupon bond with maturity T^r. 

To summarize, in the idealistic framework, we do not need to consider 9 > 
0) adapted hedging strategies, where 0 > 0) is the hltration generated by the 
observations (short rates, derivative prices, etc.) in continuous time, and we may 
restrict the admissible strategies to be adapted to the hltration generated by the 
observations at times {f, Tq, Ti, • • ■ , Tat}. 

2.2 Hypotheses on markets with liquidity costs 

We now consider markets with liquidity costs and need to precise our liquidity cost 
model. In all the sequel T_i denotes t. 

Hypothesis 2. 

We assume that, for all —1 < j < i < N, the number 7r{j,i) of zero-coupon 
bonds with maturity T* bought or sold at time Tj is measurable with respect to 
the hltration generated by {Rt, Rto^ '' iRt)- That means that the admissible 
strategies do not depend on the evolution of the rate Rg between two tenor dates 
Tm and 

Denote by T(T, U, tt) the buy or sell price for tt zero coupon bonds. In perfectly 
liquid markets, T(T, U, tt) is the linear function B{T, U)7r, where B{T, U) is dehned 
in Dehnition 1. In the presence of liquidity costs, T(T, U, tt) becomes a non-linear 
function of tt. 

Hypothesis 3. 

For all T and U, the price T(T, f/, tt) is a increasing, convex one-to-one 

map of TT from M to M, and T(T, f/, 0) = 0. 

Under the preceding hypothesis, the function T is positive when tt > 0 and 
negative when tt < 0. 

In the context of the swap we set 

T„(7r);=T(T„T,,7r) (4) 

and we only consider self-hnancing strategies, that is, satisfying 

V0<j<iV-l, Tr(k,j) + P{j) = Y (5) 


4 



2.3 Optimization objective 

In the presence of liquidity risk, the market is not complete any more and the 
practitioners need to build a strategy which minimizes a given function S (e.g. a 
risk measure) of the hedging error. Such strategies are usually obtained by solving 
stochastic control problems. These problems require high complexity numerical 
algorithms which are too slow to be used in practice. We here propose an efficient 
and original numerical method to compute approximate optimal strategies. As 
the perfect hedging leads to a null portfolio at time T^r, we have to solve the 
optimization problem 

inf E|S(H'’')], (6) 

TTGn 

where is the terminal wealth (at time T/v) given the strategy tt in the set 11 
of admissible strategies. 


3 Hedging error minimization method in a Gaus¬ 
sian framework 

The methodology we introduce in this section is based on the two following key 
observations: 

(1) We consider strategies and portfolios with finite second moment, and thus 

optimize within for some probability measure p. The Gram-Schmidt 

procedure provides countable orthogonal bases & of the separable Hilbert 
space Our set H of admissible strategies is obtained by truncating 

of a given basis, which reduces the a priori inhnite dimensional optimiza¬ 
tion problem (6) to a hnite dimensional parametric optimization problem of 
the type inGg© E\k(0, X), where 0 is a subset of MP, X is a given random 
variable, T is a convex function of 6. 

(2) The Robbins-Monro algorithm and its Chen extension are stochastic alterna¬ 

tives to Newton’s method to numerically solve such optimization problems. 
These algorithms do not require to compute X). They are based on 

sequences of the type 


d 

6*^+1 = 6*^ — p^+i—X..y+i), (7) 

where (p 7,7 > 1) is a decreasing sequence and (X.y ,7 > 1) is an i.i.d. se¬ 
quence of random variables distributed as X. 
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We here consider the case of swap in the context of a Gaussian yield curve. 
This assumption is restrictive from a mathematical point of view but is satished by 
widely used interest rate models such as Vasicek model, Gaussian affine models or 
HJM (Heath, Jarrow and Morton) models with deterministic volatilities (see e.g. 
Gibson et al. (2010) and Musiela and Rutkowski (2005)). In Babbs and Nowman 
(1999) it is shown that using a three dimensional Gaussian model is sufficient to 
£t the term structure of interest rate products. 


3.1 Step 1: finite dimensional projections of the admissible 
controls space 


Gonsider a Gaussian short rate model {Re, 6 > 0) : we either suppose that the 
dynamics of the short rate model is given, or that it is deduced from a forward rate 
model such as in the HJM approach for term structures: see e.g. Eq (6.9) in Gibson 
et al. (2010) under the additionnal assumption that the forward rate volatilities 
are deterministic. In all cases, the resulting bond price model is log-normal. 

In view of Hypothesis 2, each control 7r{j, i) belongs to the Gaussian space gen¬ 
erated by {Rt, Rtq, ■ ■ ■ , Rtj) or, equivalently, to a space generated by i{j) +1 stan¬ 
dard independent Gaussian random variables (with £{j) = j 

for one-factor models, i{j) = 2j + 1 for two-factor models, etc.) and Rt- An 
explicit orthonormal basis of the space generated by G^^\ ■ ■ ■ , is 



,/n^. 


(no,- ,n^(j))eN^O')+l 


where {Hn, n > 0) are the Hermite polynomials 


( 8 ) 


Hn{x) 


„2 /O (I I 


(_l)ngxV2^g-^V2 


dx^‘ 


(see e.g. Malliavin, 1995, p.236). 

Thus, the quantities of zero-coupon bonds bought by the trader can be written 


as 


hi) 


51 a”“(j.*)n 




m=0 


■y/UrG 


’ ^(i) ■ ■ ■ ) '^i(j))i 


(9) 


where the inhnite sum has an limit sense. A strategy can now be dehned as a 
sequence of real numbers a'^(i^{j,i) for all — 1 < j < i — 1 and n^) G 
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In order to be in a position to solve a finite dimensional optimization problem, 
we truncate the sequence [j, i)). Then a strategy is defined by a finite number 
of real parameters i), —1 < j < i < iV — 1} where, for all j, tiq) belongs 

to a finite subset of xhe truncated quantities of zero-coupon bonds 

bought by the trader write 

.U,i)= Y. n (10) 

At ^ 

We discuss the efficiency of this truncation and its convergence in Sections 3.5 
and 4.1. 

To simplify, we denote by ck = in parameters to optimize in 

(where the dimension p is known for each truncation = —1, 0, ■ • • , N — 

1)), by 77 (ck) (or, when no confusion is possible, simply tt) the hedging strategy 
corresponding to a vector a G see (10). Given the strategy 77 = 77 ( ck ), the 
terminal wealth (at time Tn) satisfies 

= E n{cx){j,N) + P{N)= '^U,N) + P{N). (11) 

— l<j<N — 

The problem (6) is now formulated as: find (x* in MP such that 

E[S («'<“•')] = inf E [S (iy>“>)] ■ (12) 


3.2 Step 2: stochastic optimization 

Using the self-financing equation (5) one can express as a function of ot and 
(( 7 ( 0 )^... , Therefore one needs to minimize the expectation of a determin¬ 

istic function of the parameter ot in and the random vector (G^°\ • • • , 

Such problems can be solved numerically by classical stochastic optimization al¬ 
gorithms, such as those introduced in the pioneering work of Robbins and Monro 
(1951) and its extensions (e.g. Chen and Zhu, 1986). We refer the interested reader 
to the classical references Dufio (1997); Chen (2002); Kushner and Yin (2003). 

In our context (12), the Robbins-Monro algorithm (7) works as follows. Start 
with an arbitrary initial condition ckq in MP. At step 7 -l- 1, given the current 
approximation cx.^ of the optimal value a*, simulate independent Gaussian ran¬ 
dom variables • • • , and compute the terminal wealth Then, 

update the parameter ct by the induction formula 

= a., - p^+i [5(IU(“-))] , (13) 
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where (p^) is a deterministic decreasing sequence. In addition, one can use an 
improvement of this algorithm due to Chen and Zhu (1986). Let {Ki, / > 0) be an 
increasing sequence of compact sets such that 

Ki C Int(iLi+i) and limiL; = W, (14) 

where Int(iL/+i) denotes the interior of the set iCz+i. The initial condition ckq is 
assumed to be in Kq and we set Z( 0 ) = 0 . At each step 7 in (13), if q:.y+i G 
we set /(q + 1 ) = /(q) and go to step 7 + 1. Otherwise, that is if ^ we 

set OL^^i = ckq and Z(7 + l) = /(7) + l. This modihcation avoids that the stochastic 
algorithm may blow up during the first steps and, from a theoretical point of view, 
allows to prove its convergence under weaker assumptions than required for the 
standard Robbins-Monro method. 


3.3 Summary of the method 

Our setting 


• The interest rate model satisfies: for all 0 < j < A^, there exist an integer 
^{j) and a function such that 

(Ri, Rto, • ■ ■ , Rt,) = ^j{Ru ■ ■ ■ , GW^))). 


• For all 0 < j < A^, a hnite truncation set C is given and A*^ = 

{ 0 }. 


• A strategy tt = 77 ( 0 ) is defined by ck = (a^'O)(j,z), — 1 < j < i < N — 
1, nQ) G A*^-^)). More precisely, the number of zero-coupon bonds with matu¬ 
rity Ti bought or sold at time Tj is 


7r(jA) 


hi) 


Ti(j)eAO) 


m=0 


g„„(GW) 


for z < A^ — 1, 


(15) 


and 71 {j, N) is deduced from the self-hnancing equation 


^7r(z,j) + P(j) = ^Tj,i(7r(j,i)) (16) 

i<j i>j 


(one possibly needs to use a classical iterative procedure to solve this equation 
numerically). 


• One is given an increasing sequence of compact sets {Ki,l > 0) satisfying 
(14) and a sequence of parameters (^^,7 > 1) decreasing to 0. 
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Our stochastic optimization algorithm 

Assume that the parameter cx^ = (j, i), —1 < j < i < iV — 1, n(j) e and 

are given at step 7. At step 7 + 1; 


1 . 

2 . 


Simulate a Gaussian vector (G® • • ■ , 

Deduce the quantities of zero-coupon bonds from (15) 

_ Tj (\ 

7i^+i{j,i)= n , 

nyjSAO) «^=0 

and get 7r^+i(j, iV) from the self-hnancing equation 


and (16); 


for i < iV — 1, 


+ P{j) = ^%(7r^+i(jA))- 

i<j i>j 


3. Compute the terminal wealth 

wAi’ = E’^-'+‘W'7') + -pw^ 

j<N 


4. Update the parameters 


^ 7+1 


= CK, 


P7+l'^ Ot 




5. If CK^+i ^ A'q, set cx^+i = ckq and /(q + 1) = /(q) + 1. 

6. Go to 1 (or, in practice, stop after V steps). 


3.4 Error analysis 


In this subsection we study the convergence (Theorem 4) and convergence rate 
(Theorem 6) of the stochastic algorithm used in Step 2, when the total number of 
steps T tends to inhnity. We introduce some notation. Recall (13) and write 


O^yT-l |^Vq;S*( )J 

Here, is given by 

= v„s(Tq“;') - E [v„s(t 77>)' ■ 


(17) 


(18) 
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The last term in (17) represents the reinitialization of the algorithm if 

^ -^/( 7 ) he. is fixed such that Q:..y+i = ckq- 

Let us now recall the convergence theorem obtained by belong (2008, Theo¬ 
rem 1 ) in our setting. 

Theorem 4. Assume 

(Al) The function cx i—>■ E[S'(1T^“^)] is strictly concave or convex, 

(^ 2 ) E7P7 = 00^ E7P7 < OO; 

(A3) The function a 1 —)■ E[|| Vct*S'(lT*'“^)|p] is bounded on compact sets. 

Then the sequence ( 0 : 7,7 > 1) converges a.s. to the unique 
optimal parameter a* such that 

inf^E = E . 

Hypothesis 3 and (15) imply that (A3) is satished. Before giving examples of 
situations where (Al) is fulhlled, let us check that IT^") is a concave function of 

o. 

Proposition 5. The terminal wealth is a concave function of the parameter 

OL. 

Proof. Recall that the terminal wealth is given by (11). The payoff of the swap 
P{N) does not depend on o. We only have to deal with the quantities 7r(z, N) 
of the zero-coupon bonds with maturity T/v bought at time T,. They satisfy the 
self-hnancing equation (5) and thus 

7r{i,N) = ^ (7r(q j)) + P(i) j , (19) 

\j<i i<j<N / 

where d'-is the inverse of the price function Tj^Ar (see (4)). Moreover, quantities 
7i{i,j), i < j < N — 1 are linear in cx (see (9)). 

Recall that Tjj is convex, thus —Tjj is concave and the argument in (19) is a 
concave function of cx. Finally is an increasing concave function, from which 
n{i,N) is a concave function of cx. □ 
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The preceding observation shows that (Al) is satisfied when S' is a utility 
function (and thus increasing and concave) and satisfies S(0) = 0. Notice that the 
optimization problem (12) then penalizes the losses and promotes the gains. In 
Sections 4 and 5 we will see another situation where Theorem 4 applies. 

Given suitable functions S, Theorem 4 guarantees the convergence of our algo¬ 
rithm towards the optimal parameters. The following theorem provides the rate 
of convergence (belong, 2013). 


Theorem 6. Let 


(v2 + 'y)^' 


( 20 ) 


for some positive Vi, V 2 and (3 G (1/2,1). Denote by the normalized centered 
error 


= 


ot^y — a* 
\//A 


Assume 


(Al) The function a is concave or convex. 

(Af) For any q > 0, the series 

^ ^ P'Y+ldM'Y+l'^{\oc^-a*\<q} 

7 

converges almost surely. 

(A5) There exist two real numbers Ai > 0 and A 2 > 0 such that 

supE < 00 . 

7 

(A 6) There exists a symmetric positive definite matrix S such that 

(A7) There exists p > 0 such that Mn > 0, d{oL*, dKn) > p, 

where dKn denotes the boundary of Kn. Then, the sequence (A.y ,7 > 1) converges 
in distribution to a normal random variable with mean 0 and covariance matrix 
linearly depending on S. 

Remark 7. As explained in detail in belong (2013, Sec. 2.4), the assumptions of 
Theorem 6 are satisfied as soon as 
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• There exists ^3 > 0 such that 


VC > 0, 


E 


sup 

l«l<c 


< CX). 


• The function ex. 1—)■ E[S'(hh^“)] is strictly concave or convex. 

The first condition is usually easily satished, e.g. when S has a polynomial growth 
at inhnity because the moments of the wealth process are typically hnite. The 
both properties are e.g. fulhlled by the example studied in Sections 4 and 5 (see 
the discussions at the beginning of these sections). 

3.5 Performance of the optimal truncated strategy without 
liquidity cost 

The numerical error on the optimal wealth decreases when the ’s tend to 
In this subsection, we provide a theoretical estimate on the error resulting from the 
truncation in (9) in the idealistic context of no liquidity cost and general Gaussian 
affine models (Dai and Singleton, 2000). 

In Dai and Singleton (2000), general Gaussian affine models are introduced 
for which, for any two times s < t, there exist standard independent Gaussian 
random variables and real numbers p, Aq, ■ ■ ■ , Xm such that the 

prices of zero-coupon bonds have the form 

B{s, t) = exp (-P - -AmG'(^)) . 

A control of the error of truncation is given in the following proposition. 

Proposition 8. In the above context, if the truncation set defined in (10) is : = 
{u-o ■ -1- < d}, then 


fd+l 

where Co and Ci are some positive constants. 

The proposition is a straightforward consequence of (11) and the next lemma 
applied to X = 7r{cx*). This lemma also allows one to precise the values of Co and 
C'l. 
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Lemma 9. Consider the random variable 


( M 

m=0 

where /i, Aq, • • ■ , Am cife real numbers, and , G^^^ are independent stan¬ 

dard Gaussian random variables. Consider the projection of X on the subspace 
of L‘^{G^^\ • • ■ , generated by 




HnJG(^'>) 

y/nj. 


no H-h nM 



We have 


||X — X '^||2 ^ exp (/i + Aq + ■ ■ ■ A^) 


d-\-l 

(Aq H-h Xm) 2 

\/('^ + 1 )! 


( 22 ) 


We postpone the proof of this lemma to the Appendix. 


4 Numerical validation of the optimization proce¬ 
dure: an example without liquidity cost 

In this section we stndy the accnracy of onr algorithm in the no liquidity cost case 
where a perfect replication strategy is known (see Section 2.1). We minimize the 
quadratic risk measure of the hedging error. 

The bond market model is the Vasicek model which is the simplest Gaussian 
model; 

dRe = A(roo — Re)d6 + adBg, (23) 

where A is the mean reverting rate, r^o is the mean of the equilibrium measure, a 
is the volatility and {Bq, 0 > 0) is a one-dimensional Brownian motion. 

Notice that 

pV 

Wu<v, = Too + (i?« - -1-/ e^^^-'^'>dBe. (24) 

J U 

Therefore, there exists an i.i.d. sequence {G^^\ G^^\ • ■ ■ , of A/'(0,1) Gaussian 
random variables such that 

VA; = 0,-- - ,iV, = $(Tfc-Tfc_i,i?T,_,,G('=)), (25) 
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where 


$(?7, r, G) = Too + e - r^) + Ga 

In our numerical experiments, we have chosen the following typical values of the 
parameters A = 10%, Too = 5%, a = 5%. With this choice of parameters, the 
mean yearly interest zero-coupon rates with maturity less than 10 years take values 
between 3% and 5%. 

Our numerical study concerns the minimization of the quadratic mean hedging 
error which corresponds to the choice S{x) = in (12). This choice penalizes 
gains and losses in a symmetric way and aims to construct a strategy as close as 
possible to the exact replication strategy. 

In the no liquidity cost case, the terminal wealth is a linear function of the 
parameter cx and therefore assumption (Al) of Theorem 4 is obviously satished. 

Given a degree of truncation d, the set is chosen as 

Jsff ;= (nQ) = (no, ••■,%)€ no H- nj < d}. (26) 

We have to optimize the real-valued parameters (j, i) for j < i and n^) G 

The quantities of zero-coupon bonds to exchange are given by (10). The choice 
of the sequence (p.y ,7 > 1) in (13) is crucial. Choose as in (20). We discuss 
the sensitivity of the method to the parameters ui, V 2 , /3 in Section 4.2.2. We also 
discuss the sensitivity of the results to the number T of steps. 

In all the sequel, we use the following notation. 

Notation For all vector cx = (o^'O) (j, i); — 1 < j < i < — 1), we set 

n((a) :=E[IF(“)]^ (27) 

where the expectation is computed only with respect to the Gaussian distribution 

4.1 Empirical study of the truncation errors (Step 1) 

In this subsection we develop an empirical validation of the projection step pre¬ 
sented in Section 3.1 

We observe that the quadratic mean hedging error decreases very fast to 0 
when the degree of truncation increases. For a notional equal to 1, the error is of 
the order of one basis point (a hundredth of percent) for a degree d = 3 and a 
small number of dates N, and for d = 4 and for larger values of N. 
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Figure 1 below shows u(q:*’'^), with the optimal parameter corresponding 
to the truncation set (26). We have used the explicitly known hnite dimensional 
projections of the optimal strategies without liquidity cost to obtain and a 
Monte Carlo procedure to compute v. Table 1 shows some values used to plot 
Figure 1. 
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0.0001 

1e-06 

1e-08 

1 e-10 

1e-12 
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1e-16 

1e-18 

1e-20 
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N 



Figure 1 


4.2 Empirical study of the optimization step (Step 2) 

The stochastic algorithm converges almost surely to the optimal coefficient a*. In 
this part, we empirically study the convergence rate in terms of the number of 
steps F and the choice of the sequence (p-y). 
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degree d 

N=2 

N=3 

0 

5.2 E-6 

3.0 E-5 

1 

5.4 E-9 

3.1 E-8 

2 

3.7 E-12 

2.0 E-11 

3 

1.9 E-15 

1.9 E-14 

4 

2.2 E-18 

3.9 E-15 


Table 1: 


4.2.1 A typical evolution of (o^) 

In this subsection, we consider a swap with two payment dates {N = 2). We 
consider the truncation set = {0,1}. The objective is to approximate a* = 
(q;°’*(—1, 0), 1,1), q;°’*( 0, 1), Q;^’*(0,1)). In Figure 2, the four parameters a = 

(q;°(—1, 0), q;°(—1,1), q;‘^(0, 1), q;^( 0, 1)) evolve according to (13) where the sequence 
(^ 7,7 > 1) is dehned by (20) with vi = 10’^, V 2 = 1 and 13 = 1. 

As expected, the sequence (0:7) converges to ck*. However, the evolution is 
quite slow although we have empirically chosen the parameters vi, V 2 and /3 in a 
favorable way. 

In Figure 3, we plot (in violet) n(Q;°(—1,0), 1), q;°’*(0, 1), q;^’*( 0, 1)) as a 

function of q;°(— 1,0) and q;°(— 1,1). We also plot in green the path (n(Q:.y),0 < 
7 < F). The hgure shows that after F = 10000 steps the hedging error is small 
though the optimal parameters have not been approximated accurately (notice 
that the violet surface is flat). 
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Figure 2: Evolution of the parameters a°(—1,0), a°(—1,1), a°(0,1), q;;^( 0, 1) in 
terms of 7. 
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-0.03 

-0.035 
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-0.005 

- 0.01 

-0.015 

- 0.02 

-0.025 

-0.03 

-0.035 


Figure 3: Violet surface: —u(.,a°’*(0,1), 1)). Green curve: evolution 

of (q;° 7 (— 1 , 0 ), q;°(—1 ,1 ), —u(q;°(—1 , 0), Q;°(— 1 ,1 ), a°’*(0,1), 1)) in terms of 

7 = 1 , • • • , 10000 (ui = 10 , U 2 = 1 and /? = 0 . 6 ). 
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4.2.2 Sensitivity to the choice of the sequence (p^) 

Theorem 4 states the convergence of the optimization method for all sequence (p^) 
satisfying (A2). We here study the sensitivity of the results to the parameters ui, 
^ 2 , /d of sequences of type (20) and to the total number of steps T. 

Tables 2 and 3 show the expected value function obtained after T = 10i?4, 
10i?5 and 10i?6 steps. The expected value function is estimated by means of a 
classical Monte Carlo procedure. Table 4 shows the same results with a sequence 
= Vi which does not satisfy condition (A2). 


r 

1 

10 

100 

1000 

10 000 

20 000 

10E5 

10E6 

10E4 

7.2 E-4 

2.6 E-6 

1.1 E-6 

8.9 E-7 

5.1 E-4 

8.3 E-4 

1.0 E-1 

1.3 E-1 

10E5 

1.4 E-5 

9.3 E-7 

8.9 E-7 

4.7 E-7 

2.3 E-8 

1.7 E-2 

9.1 E+3 

9.5 E+5 

10E6 

3.8 E-6 

9.2 E-7 

7.7 E-7 

1.5 E-7 

6.6 E-12 

5.4 E-12 

1.0 E-5 

15.8 


Table 2: ^(cKr) (^2 = 1000, /3 = 0.6) 


r 

1 

10 

100 

1000 

10 000 

13 000 

2E4 

10E4 

8.8 E-3 

1.0 E-3 

7.8 E-6 

9.8 E-7 

8.8 E-7 

1.7 E-6 

9.7 E-7 

10E5 

6.6 E-3 

1.1 E-4 

9.3 E-7 

9.0 E-7 

6.8 E-7 

5.5 E-7 

4.8 E-7 

10E6 

4.5 E-3 

1.8 E-5 

9.3 E-7 

8.7 E-7 

5.1 E-7 

3.6 E-7 

2.9 E-7 


r 

1E5 

5E5 

1E6 

2E6 

3E6 

4E6 

5E6 

10E4 

1.3 E-3 

3.3 E-2 

7.3 E-1 

4.5 E+1 

1.1 E-2 

6.5 E+2 

6.9 E+1 

10E5 

6.6 E-7 

1.2 E-5 

4.1 E-4 

2.7 E-1 

8.1 E-3 

1.7 E-2 

1.4 E+1 

10E6 

6.8 E-8 

1.1 E-10 

5.3 E-12 

7.4 E-12 

7.0 E-12 

2.6 E-5 

1.5 E-6 


Table 3: ^(cKr) (^2 = 1000, /3 = 0.9) 


Vl 

r 

1 

2 

4 

6 

8 

10 

12 

20 

10E4 

7.6 E-7 

6.9 E-7 

2.5 E-7 

2.9 E-7 

1.4 E-6 

1.9 E-7 

7.6 E-7 

6.1 E-6 

10E5 

6.9 E-7 

6.8 E-7 

4.1 E-7 

1.3 E-7 

2.9 E-7 

3.2 E-7 

3.9 E-6 

3.2 E-4 

10E6 

3.0 E-8 

1.0 E-9 

5.1 E-12 

4.3 E-12 

5.1 E-12 

4.2 E-12 

5.9 E-12 

7.8 E-6 


Table 4: ^(cKr) for a constant sequence = Vi 

We observe that the efficiency of the algorithm depends on the choice of the 
parameters ui, ^ 2 , and is really sensitive to it when the total number of steps T 

is small. 
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When r becomes large (e.g F = 10ii^6), then the algorithm may seem to diverge 
if /3 is chosen carelessly. In fact, as the sequence (p.y ,7 > 1) satisfies hypothesis 
(A2) of Theorem 4, the algorithm converges to the optimal parameters but it is far 
from (X* after 10i?6 steps. However, for each value some vi reduces the mean 
square hedging error to 5 E-12. 

5 An empirical study of the bid-ask spread costs 
impact 

We here present numerical results corresponding to two piecewise linear liquidity 
cost functions T; 



(28) 


(29) 


Despite the fact that we know there is no perfect hedging strategy in this 
context, we suppose the holder receives a null cash at time t (which is the price of 
the swap in a no liquidity cost market). 

We now shortly check the convergence of the algorithm. Given piecewise linear 
cost functions T, it is easy to prove that the terminal wealth is piecewise 

linear in cx (see the proof of proposition 5). Therefore, assumption (Al) of The¬ 
orem 4. However, Theorem 4 does not apply to our context since 4/)^ and 
are piecewise linear and therefore are not continuously differentiable everywhere. 
Replace and ft'l by smooth approximations obtained by convolutions with ker¬ 
nels of the type 1/\/2T^e exp(—(2£)), e small. Let cx* be the unique optimal 
parameter corresponding to the new cost functions (existence and uniqueness of 
CK* are provided by Theorem 4). In view of Rockafellar and Wets (1998, Th 7.33) 
(cK*) tends to cx* when e tends to 0. 

The preceding consideration to prove convergence is more theoretical than prac¬ 
tical: in practice, the numerical results do not differ when e is small or e is null. 

Finaly, notice that the piecewise linearity in cx of implies that conditions 
in Remark 7 are satisfied, so that Theorem 6 precises the convergence rate. 
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5.1 Taking liquidity costs into account is really necessary 

Consider two different strategies: (i) the strategy corresponding to the optimal 
parameters ckq in the idealistic model withont liqnidity costs and (ii) the null 
strategy <5o defined as 

<^0 (j, 0=0) for all —l<j<i<N — 1 and e 

To satisfy the self-hnancing assumption (5), at time Tj the payoff P{j) of the swap 
(2) is used to buy zero-coupon bonds with maturity T^. 

Figure 4 shows —^(cko) and —n(5o) in terms of the parameter A where the cost 
function 4/ is as in (28). The mean square hedging error dramatically increases 
when, in the presence of liquidity costs, the trader uses the strategy which is 
optimal in the no liquidity cost context. When the liquidity cost A is larger than 
4%, it is even worse to use this strategy than to use the 5o strategy! 


0.02 



1 


1 

1 


1 

o o 

Q to 

'> > 

) . 

/ 

/ 

/ 

1 _ 

/. 

/ 

/ 

/ 

/ 










\ 

\ 

\ 

\ 

\ 









.\" 

N 

\ 

\ 

\ 









\ 

\ 

\ 

. 








\ 


_1 

_ 

_1 

\ 

\ 

\ 

_ 

_ 

_1 

_1 

_ 


- 0.02 - 


- 0.04 - 


- 0.06 - 


- 0.08 - 


- 0.1 


0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 


Figure 4: (short dashes) —n(5o) (long dashes) —^(cko) in terms of the liquidity 
cost A 
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5.2 Probability distribution of the hedging error in the case (28) 

In this section, the liqnidity cost fnnction d' is chosen as in (28). 

After r steps of the stochastic optimization procednre with a sample uj of the 
Gaussian vector ■ ■ ■ , 7 = 1, • ■ ■ , F), one obtains a random approxi¬ 

mation Q:r(i^) of the optimal parameter o*. 

Figure 5 shows the probability distribution of the random variable n(Q:r(i^)) 
for r = 10000 and Table 5 shows its mean and standard deviation for different 
values of A. 


A 

0 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

Mean 

2.5E-32 

0.0031 

0.012 

0.024 

0.038 

0.049 

0.058 

0.065 

0.11 

0.12 

Std dev. 

6.7E-33 

4.2E-5 

2.8E-4 

1.7E-3 

0.016 

0.017 

0.046 

0.081 

1.1 

1.3 


Table 5; Empirical mean and standard deviation of ^(ckioooo)- 
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Density . 
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-0.1 -0.08 -0.06 -0.04 -0.02 0 


Figure 5: Empirical density of —n(Q:r(ti;)) (ui = ^2 = /3 = 1, F = 10000, A = 0.04). 
The vertical line corresponds to the v{a*). 


5.3 Hedging error in the case (29) 

In this section, the liquidity cost function T is chosen as in (29). 

Figure 6 shows —^(cKr) for F = 0 (green) and F = 10E6 (violet). The initial 
parameter ckq is the optimal one for a market without liquidity cost. 
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We observe that the main part of the loss is saved thanks to the optimization 
procedure. In Figure 7, we zoom on the surface resulting from the optimization 
procedure. Notice that the cost function (28) is equal to the cost function (29) 
in the particular case (7 = 0. So, Figure 7 allows to compare the value functions 
corresponding to these two cost functions. 



Figure 6; Violet surface: —u(Q:r) (F = 10F^6, vi = V 2 = 100, /I = 0.6) for a 
liquidity cost function (29) in terms of the value A of the liquidity cost and the 
size C of the compact. Green surface:— u(q:o). 


5.4 Influence of the initial value ckq of the optimization pro¬ 
cedure 

In Figure 8, we draw two functions of the liquidity cost A: —^(cKr) and —u(5r), 
where (or) and (5r) are the parameters obtained after F steps of the optimization 
procedure but with different initial values ckq and <5o as described in Subsec. 5.1. 
The performance of the strategies obtained after F = 10F^6 steps are quite similar. 
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Figure 7; A zoom on the violet surface of Figure 6 


It means that the sensitivity to the arbitrary initial parameter ckq is not observable 
any more after F = 10776 steps. 

5.5 Reducing the set of admissible strategies 

Recall Hypothesis 2. So far, our admissible strategies at time Tj depend on all 
the past and present rates Rt-, - ■ ■ , Rry Thus the number of parameters [j, i) 
to optimize is at least of the order of magnitude of the binomial coefficient 
where N is the number of dates and the degree of truncation d is dehned as in 
(26). This order of magnitude is a drastically increasing function of N. This 
crucial drawback leads us to try to simplify the complexity of the control problem 
(6) by reducing the size of the set of the admissible strategies H. Observe that 
the optimal strategy under the perfect liquidity assumption has the property that 
71* (j, i) only depends on Rt^ . This observation suggests to face large numbers of 
dates by reducing the set of controls to controls depending only on a small number 
of recent interest rates Rt^Rt^.i, ■ ■ ■ , Rtj^^- 
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Figure 8: (short dashes) —u(<5r) (long dashes) —^(cKr), in terms of the liquidity 
cost A (r = 10F^6, Vi = 0.1, V 2 = 100, /3 = 0.6, d = 3). 


Figures 9 and 10 illustrate that the optimal control problem ( 6 ) with admissible 

strategies defined as in hypothesis 2 may be used as benchmarks to solve control 
problems. Consider swaps with N = 5 and iV = 10 dates of payment. In each 
one of these two cases, we study the effect of choosing q = 0 (that is at time Tj, 
admissible strategies only depend on Rt^), <1 = ^ (admissible strategies depend on 
Rt^ and i?T,_i), g = 2 . 

Figure 9 shows the performance of the corresponding strategies obtained after 
F = 10i?6 steps of the optimization algorithm for a swap with N = 5 dates of 
payments. 

Figure 10 shows similar quantities for a swap with iV = 10 dates of payment. 
We observe that the numerical computation of the optimal strategy is quite un¬ 
stable when q is too big, which reflects the difficulty to solve a high dimensional 
optimization problem. Therefore, one necessarily must choose q small in order 
to get accurate approximations of optimal strategies belonging to reduced sets of 
admissible strategies. 
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Figure 9: —^(cKr) in terms of A for strategies depending only on recent rates 
Rt^, • • • , RTj.q {q = 0,1, 2, r = 10i?6, vi = 0.1, V 2 = 100, jS = 0.6, d = 3) for a 
swap with N = 5 dates of payment. 
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Figure 10: —^(cKr) in terms of A for strategies depending only on recent rates 
Rt^, • • • , RTj.q {q = 0,1, 2, r = 10i?6, vi = 0.1, V 2 = 100, (3 = 0.6, d = 3) for a 
swap with = 10 dates of payment. 
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6 Conclusion 


Stochastic control problems generally have no explicit solutions and are difficult 
to solve numerically. In this paper, we have proposed an efficient algorithm to 
approximate optimal allocation strategies to hedge interest rate derivatives subject 
to liquidity costs. 

As discussed above, our methodology is constructive and efficient in a Gaussian 
paradigm. We project the admissible allocation strategies to the space generated 
by the hrst Hermite polynomials and use a classical stochastic algorithm to op¬ 
timally choose the coefficients of the projection in order to optimize an expected 
function of the terminal hedging error. 

We have illustrated this general approach by studying swaps in the presence of 
liquidity costs. We have discussed the performances of the numerical method in 
terms of all its algorithmic components. 

We emphasize that our methodology can be applied to many control problems, 
e.g., the computation of indifference prices, when the model under consideration 
belongs to a Gaussian space. 
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7 Appendix 


Proof of Lemma 9 

Proof. Bet us hrst prove the result for M = 0, that is X = exp [p + XG). 


X 


+00 

^a(n) 

n=0 


Hn{G) 

\/n\ 
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where 


a{n) = E 


X 


Hn{G) 

\/n\ 


= exp(/i)E 


exp(AG) 


Hn{G) 

y/n). 


We then use the identity: a /2 _ ^Hj {x) and obtain 

a{n) = exp(/i + AV2)E | ^ 

, J- vn\ 


As {Hj{G)/\/^.,j > 0) is an orthonornial basis of L‘^{G) we have 


a{n) = exp(/i + A^/2) 


A^ 


y/n^. 


Let be the projection of X on the subspace generated by the Hq, ■ ■ ■ , Hd : 

TT / d 


^ TT (/^\ ^ \n TT //^\ 

X^ = y a{n)ytl = V 


X - X'^ = y 

The truncation error is 


A- Hr,{G) 


n=d-\-l 


n\ vn\ 


^ \2n 

\\X - X% = y ^ 


n=d+l 


< e 


(d + iyX^ „! 


2n 


< e 


2fi+2X 


2 


(d + iy: 


The desired result thus holds true for M = 0. Let M be a positive integer. We 
have 


M 


iix - x'^wi = e2^+^o+'"+^" y JJ 


A 


2 ?T- j 7 


Recall the classical identity 


(oq + ■ ■ ■ + O-l) 

xT 


N 


«-oH- hnM>d m=0 


E n 


[Ur. 


noH- \-n]^=N m=0 


.nr. 


(30) 
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Thus, 


M 


E n 

noH- hnM>d m=0 


X 


‘2iTtr) 


M 


A Ml 


E E n 

N=d+1 noH- \-nM=N' m=0 


X 


2nn 


irir. 


E 


(Aq H-h Xm) 


N 


N\ 

I 2 'i'i+1 


N=d+1 

< + + exp (A^ + 


(d+1)! 


This ends the proof for all positive M. 


-^ -^m) • 
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